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Neutrinos offer a particularly promising eye on the extreme Universe. Neutrinos 
are not attenuated by intervening radiation fields such as the Cosmic Microwave 
Background, and so they are messengers from the very distant and very young 
phase of the universe. Also, neutrinos are not deflected by cosmic magnetic fields, 
and so they should point to their sources. In addition, there are particle physics 
aspects of neutrinos which can be tested only with cosmic neutrino beams. After a 
brief overview of highest-energy cosmic ray data, and the present and proposed ex- 
periments which will perform neutrino astronomy, we discuss two particle physics 
aspects of neutrinos. They are possible long-lifetime decay of the neutrino, and a 
measurement of the neutrino-nucleon cross-section at a CMS energy orders of mag- 
nitude beyond what can be achieved with terrestrial accelerators. Measurement 
of an anomalously large neutrino cross-section would indicate new physics (e.g. 
low string-scale, extra dimensions, precocious unification), while a smaller than 
expected cross-section would reveal an aspect of QCD evolution. We then discuss 
aspects of neutrino-primary models for the extreme-energy (EE) cosmic ray data. 
Primary neutrinos in extant data are motivated by the directional clustering at 
EE reported by the AG AS A experiment. We discuss the impact of the strongly- 
interacting neutrino hypothesis on lower-energy physics via dispersion relations, 
the statistical significance of AGASA directional clustering, and the possible rele- 
vance of the Z-burst mechanism for existing EE cosmic ray data. 

1 Introduction 

Detection of ultrahigh-energy neutrinos is one of the important challenges of the 
next generation of cosmic ray detectors. Their discovery will mark the advent of 
neutrino astronomy allowing the mapping on the sky of the most energetic, and 
most distant, sources in the Universe. In addition, detection of extreme-energy (EE) 
neutrinos, those at 10 20 eV and beyond, may help resolve puzzles associated with 
the giant air-shower events observed with energies beyond the Greisen-Zatsepin- 
Kuzmin (GZK) limit of Eqzk = 5 x 10 19 eV. In this context, neutrino observations 
may validate Z-bursts, topological defects, superheavy relic particles, new strong- 
interactions, etc. 

Measurement of the neutrino cross-section itself at EE is of considerable impor- 
tance to particle physics. Our highest-energy knowledge of the neutrino-nucleon 
cross-section comes from the HERA experiment, at a CMS energy of 0.2 TeV. 
The CMS energy for a CR of energy E = 10 20 E20 eV on an air nucleus is 
y/s = 0.5 A/E20 PeV. This can be compared to the values of terrestrial hadron 
accelerators, 2 TeV at Fermilab, and 14 TeV to occur at CERN's LHC. Further- 
more, if the primary CR is a neutrino, its energy is not shared among partons, so 
its "reach" is even larger than the comparison with Fermilab and the LHC indi- 
cates. Estimates of the cross-section at 10 20 eV require QCD extrapolations over 
three orders of magnitude beyond HERA in CMS energy. Thus, a measurement of 
the cross-section at 10 20 eV will test QCD with a large lever arm. Additionally, 
there may well be new non-SM physics revealing itself in the neutrino sector at 
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EE, arising from new physics thresholds not accessible to terrestrial accelerators. 
Needless to say, the value of the cross-section is also crucial to the future evolution 
of "neutrino telescope" detectors, for the neutrino event rate is proportional to this 
cross-section (more on this later). 

And finally, we note that the stability of the neutrino is best tested with PeV 
cosmic neutrinos. The "cosmic" aspect guarantees a long decay path, while the 
PeV energy is a compromise between minimizing the boost factor 7 = E/m and 
ensuring events from extragalactic rather than atmospheric sources. 

Several sources of neutrinos with PeV to ZeV (10 21 eV) energies are possible, 
ranging conservatively from AGNs and GRBs, to exotic top-down production. The 
latter may in principle even provide energies up to the grand-unified mass of ~ 
10 24 eV. A nice review of sources, classified according to their speculative nature, 
was given a few years ago by Protheroe. El Some of the most recent ideas for EE 
neutrino sources are discussed in our later section on Z-bursts. 

In addition to the flux of neutrinos produced in cosmic engines, there is also 
a reasonably guaranteed prediction for a flux of "GZK neutrinos" in the energy 
range 10 15 to 10 20 eV, based on the observed flux of cosmic ray (CR) protons at 
and above the GZK limit. These neutrinos result from the decay of charged pions 
photo-produced by the interaction of super-GZK nucleons on the CMB background. 
This flux is expected to peak in the decade 10 17 to 10 18 eV for uniformly-distributed 
proton sources, and around 10 19 eV for "local" sources within ~ 50 Mpc of earth. 

On the experimental side, past observations of EE cosmic rays have been made 
by the Volcano Ranch, Haverah Park, AKENO, and Fly's Eye experiments. El The 
record energy, held by the Fly's Eye event of a decade ago, is 3 x 10 20 eV, or 
50 Joules. This truly macroscopic amount of energy equals that in a professional 
baseball pitcher's fastball, which consists of 10 27 nucleons. This latter comment 
shows that somehow Nature is 10 27 times more efficient at acceleration than is the 
human arm. 

Most of the catalogued EE events come from the AG AS A experiment, the suc- 
cessor to AKENO. The HiRes experiment, the successor to Fly's Eye, has an ex- 
posure similar to AGASA, but an EE event rate lower by almost an order of mag- 
nitude. The discrepancy between these two experiments is shown in Fig. ^ There 
is clearly a systematic error in the energy assignments of (at least) one of the ex- 
periments. It has been noted, and is somewhat evident in the figure, that even 
the spectral dip ("ankle") occurs at different energies in the two experiments. The 
AGASA experiment uses ground-based scintillator to measure showers, whereas 
HiRes uses optics to measure the near-UV light emitted by atmospheric N2 excited 
by the passing shower. Fortunately, both measurement techniques will be present 
in the much larger AUGER experiment under construction in the high plateau 
of Argentina. In this experiment, simultaneous measurement of the same shower 
with both techniques will allow a cross-check of systematics effects, thereby set- 
tling in the next two years the issue of AGASA vs. HiRes. De Marco et al. ™ give 
a quantitative appraisal of the AGASA vs. HiRes discrepancy, and the resolving 
potential of AUGER, and ultimately, of the "Extreme Universe Space Observatory" 
(EUSO) experiment. ^ It is worth noting that if the lower HiRes rate turns out to 
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be correct, and even if the GZK cutoff is realized, the region above Eqzk remains a 
fertile arena for new physics. The depletion or total absence of nucleon primaries in 
this energy region would constitute a background free playground for the discovery 
of new primary particles, such as neutrinos, supersymmetric particles pi magnetic 
monopolesP and any other exotic quanta. Moreover, the MFP for photons above 
10 21 rapidly increases to ~ 100 Mpc. This means that if there are sources of 
such energetic photons, they should be unmasked by the proposed larger aperture 
experiments. 



• HiRes-2 Monocular 
■ HiRes— 1 Monocular 

t AGASA 




log 10 (E) (eV) 

Figure 1. E 3 times the UHE Cosmic Ray Flux, with data from HiRes and AGASA (from T. 
Abu-Zayyad et al., HiRes Collaboration Also shown is a fit to the HiRes data of a uniform cosmic 
distribution of sources cutoff at E = 10 21 eV and subject to e + e — losses and the GZK photo-pion 
absorption mechanism. Not shown here is the Fly's Eye event at 3 X 10 20 eV. 



The EUSO experiment will mount a two meter Fresnel lens on the International 
Space Station (ISS), 400 km above the earth. The lens will look down on the 



weiler'YITP: submitted to World Scientific on February 2, 2008 3 



earth's atmosphere, and focus the near-UV N2 fluorescence emitted in air-shower 
cascades occurring in the moonless night sky. EUSO is scheduled for deployment 
in early 2009. Eventually, the free- flying multi-satellite experiment "OWL" may 
follow EUSO. In terms of US states, the Auger field of view may eventually equal 
the area of Rhode Island, while EUSO and OWL may equal Texas, or more. The 
target-mass of atmosphere available to EUSO is a teraton. For these space-based 
experiments, the 1/r 2 loss at the 400-500 km height limits their sensitivity to events 
with energy above ~ 10 19 eV, thereby providing a natural filter to select only the 
most extreme-energy CRs. 

The ICECUBE experiment is optimized to measure the neutrino flux at a PeV 
and above. It consists of strings of optical modules triggering on the Cerenkov light 
emitted by showers and muons in a gigaton of Antarctic Ice. This experiment builds 
on AMANDA, the proof-of-principle prototype experiment presently operating in 
South Polar ice. ICECUBE is funded, and is expected to be fully arrayed near 
the end of this decade. Prototypes for ocean detectors are also under construction. 
These are ANTARES and NESTOR in the Mediterranean Sea. A proposed gigaton 
array is named NEMO, possibly to be sited off of Sicily. 

The far future of EE neutrino astronomy may belong to radio frequency detec- 
tors. The idea here is that the amplitude for Cerenkov emission at a particular 
wavelength A coherently sums the net charge of the shower in the length A. Since 
the net charge is known to be about 20% of the total charge, and the total charge 
scales simply with shower energy as ~ E/GeV, it follows that the Cerenkov rate 
scales as (£'/GeV) 2 at long wavelength. Thus, radio wavelengths are tremendously 
favored at extreme-energy. The RICE experiment at the South Pole is a prototype 
for radio detection. The experiment ANITA may be the first to witness an EE cos- 
mic neutrino event. ANITA is a balloon experiment at the South Pole, configured 
to detect radio signals generated by showers originating from "earth-skimming" 
neutrinos interacting just below the horizontal (more on earth-skimmers in section 
OJ). Eventually, large radio arrays may be deployed with ICECUBE, or in the mas- 
sive salt domes which are known to populate the underground in the southeastern 
US. 

Another noteworthy experimental idea is the mountain-to-mountain "telescope" 
configured to see horizontal, neutrino-induced showers emanating from one moun- 
tain toward the lens located across the valley on the other. ^ A useful table of many 
happening and proposed EE neutrino experiments, assembled by Peter Gorham, is 
available at 

http :/ / astro . uchicago . edu/home /web / olinto / aspen/gorham_t able . htm 

In what follows we develop each of the theoretical themes noted earlier in this 
introduction. 

2 Astrophysical Test of Neutrino Decay 

Neutrinos from astrophysical sources are expected to arise dominantly from the 
decays of pions and their muon daughters, which results in initial flavor ratios 
<j) v<i : cf> u : <j) Ur of nearly 1:2:0. As the cosmic neutrinos travel over many 
oscillation lengths, they lose their phase information ("decohere") and arrive at 
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earth as an incoherent ensemble of mass eigenstates. The relative fluxes of each 
mass eigenstate are given by 4> Vi — J2 a 0t° urCC I 1 2 ' wnere U a i are elements of 
the neutrino mixing matrix. For three active neutrino species there is now strong 
evidence to suggest that and v T are maximally mixed and U e s ~ 0. Consequently, 
the neutrino mass eigenstates are produced in the ratios 4> Vl : 4> V2 : 4> Vz =1:1:1, 
and so arrive at earth in this same ratio. Due to the newly determined values of 
neutrino mixing parameters with their small observational uncertainties, this ratio 
is robust. 

It was recently showrU^that this robustness allows for sensitive tests of neutrino 
decay, as decay would alter the measured flavor ratios in a strong and distinctive 
fashion. Conversely, a measurement of ratios other than 1:1:1 can be construed 
as evidence for unstable neutrino mass. It was further shown that neutrino decay 
cannot be mimicked by either different neutrino flavor ratios at the source or other 
non-standard neutrino interactions. In this section we elaborate on these results. 

We restrict our attention to the two body decays 

V{ —tv,j+X and v\ — > Vj + X, (1) 

for which limits are too weak to eliminate the possibility of astrophysical neutrino 
decay by a factor of about 10 7 x (L/100 Mpc) x (10 TeV/E). OS Here, v, are 
neutrino mass eigenstates and X denotes a very light or massless particle, e.g. a 
singlet Majoron. Radiative two-body decay modes and three-body decays of the 
form v — > vvv need not be considered, as they are already constrained beyond 
what can be inferred from cosmic fluxes, by the absence of photons in appearance 
searches, and by bounds on anomalous Zvv couplings, respectively. 

We also assume that the decays are complete, i.e. none of the decaying state 
arrives at earth. This is reasonable because even the shortest distances are typically 
hundreds of Mpc, and the typical energies in a steeply falling spectrum are not too 
large. The assumption of complete decay means that the distance and intensity 
distributions of sources need not be considered. Finally, neutrinos and antineutrinos 
need not be distinguished because their cross sections rapidly approach each other 
above 10 TeV. 

First suppose that there are no detectable decay products, that is, the decay- 
ing neutrinos simply disappear. Such would be the case for decay to "invisible" 
daughters such as a sterile neutrino, or for decay to active daughters if the source 
spectrum falls sufficiently steeply with energy so that the flux of daughters with de- 
graded energy makes a negligible contribution to the total flux at that energy. Since 
coherence is lost, one has for the flavor fluxes (assuming completeness: L 3> Tj), 

^ a {E) — ^cf>»°;™(E)\Up s \ 2 \U as \\ (2) 

where the sum on s is over just the stable states. 

The simplest case (and the most generic expectation) is a normal hierarchy in 
which both v% and vi decay, leaving only the lightest stable eigenstate v\. In this 
case the flavor ratio is Ul x : : U^i- Neglecting [7 e3 = 0, one then has 

<t>v e ■ ■ <t>u T = cos 2 Q : i sin 2 Q : \- sin 2 ~ 6 : 1 : 1, (3) 
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where Oq is the solar neutrino mixing angle, which we set to 30°. In the case of an 
inverted hierarchy, v 3 is the lightest and hence stable state, and so 

4> Ve : ^ : 4> Vt = U 2 e3 : U% : C/ r 2 3 = : 1 : 1. (4) 

Interestingly, both cases have extreme deviations of the flavor ratio from the 
1:1:1 in the absence of decays, which provides a very useful diagnostic. Assuming 
no new physics besides decay, a ratio of 0„ e : <f> v greater than 1 suggests the 
normal hierarchy, while a ratio smaller than 1 suggests an inverted hierarchy. It is 
also interesting to note that complete decay cannot reproduce 1:1:1. One of the 
mass eigenstates does have a flavor ratio similar to 1 : 1 : 1, but it is the heavier 
of the two solar states and cannot be the lightest, stable state. (A possible but 
unnatural exception occurs if only this state decays) . 

These clear and striking predictions depend strongly on the recent progress 
which determined the neutrino mixing parameters. In particular, it is very signifi- 
cant that 6*0 ~ 30° is well below the maximal 45° , for which Eq. would instead 
be a much less dramatic 2:1:1. In addition, Oq < 45° means that 5m\ 2 > 
and hence that v 2 (with flavor ratios 0.7 : 1 : 1) can never be the lightest mass 
eigenstate. Maximal atm and very small C/ e3 also make the predictions clearer. 

Quite different ratios may result, depending on which of the mass eigenstates 
are unstable, the decay branching ratios, and the hierarchy of the neutrino mass 
eigenstates. Also, when the appearance of daughter neutrinos from the decay cannot 
be neglected, the equations are more complicated, but only slightly so. That case 
is also treated in Beacom et al,E3for both Dirac and for Majorana neutrinos. I do 
not treat that case here, but instead, list some possibilities for the normal hierarchy 
in Table I. 



Table 1. Flavor ratios for various decay scenarios (normal hierarchy only). 



Unstable 


Daughters 


Branchings 


<f>u e ■ 0i/„ : <f>v T 




anything 


irrelevant 


6:1:1 




sterile 


irrelevant 


2:1:1 


v-i 


full energy 
degraded (a = 2) 


B%^2 = 1 


1.4: 1 : 1 
1.6 : 1 : 1 


v-i 


full energy 
degraded (a = 2) 


B^i = 1 


2.8 : 1 : 1 
2.4 : 1 : 1 




anything 


B 3 ^ = 0.5 
B 3 ^ 2 = 0.5 


2:1:1 



An important issue is whether there are other scenarios (either non-standard 
astrophysics or non-standard neutrino properties) that would give similar ratios. 
The answer is that since the mixing angles 0q and 6* a t m are both large, and since 
the neutrinos are produced and detected in flavor states, no initial flavor ratio can 
result in a measured tf> Ve : (j) Ufl ratio anything like that of our two main cases, 6 : 1 
and 0:1. In terms of non-standard particle physics, decay is unique in the sense 
that it is "one-way" . Oscillations or magnetic moment transitions are, on the other 
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hand, "two-way". Since the initial flux ratio in the mass basis is 1 : 1 : 1, magnetic 
moment transitions between (Majorana) mass eigenstates cannot alter this ratio, 
due to the symmetry between i — ► j and j — ► i transitions. However, if neutrinos 
have Dirac masses, magnetic moment transitions (both diagonal and off-diagonal) 
turn active neutrinos into sterile states, so the same symmetry is not present. Yet 
the process will average out at 1/2, so there is no way to leave only a single mass 
eigenstate. 

Experimentally, the number of muon tracks and the number of showers (charged- 
and neutral-current combined) are accessible. The relative number of shower events 
to track events can be related to the most interesting quantity for testing decay 
scenarios, i.e., the v e to ratio. The precision of the upcoming experiments, e.g. 
ICECUBE, should be good enough to test the extreme flavor ratios produced by 
decays, and thereby discover or limit neutrino decay. 

3 Can't Lose Theorem for Smaller/Larger Neutrino Cross-Section 

The expected rates for neutrino observation in approved and proposed experiments 
are proportional to ovjy. The grossly-extrapolated cross section at 10 20 eV is ~ 
10 _31 cm 2 . If Nature offers a smaller cross section, then the main detection signal 
proposed for UHE neutrino experiments would be compromised. On the other hand, 
the extrapolated cross-section may be too low, for it ignores possible contributions 
from new physics that may enter in the ~TeV to PeV scale accessible to CR physics 
but inaccessible to terrestrial accelerators. 

Consider, for example, the space-based experiments EUSO and OWL. The event 
rate for nearly horizontal air showers (HAS) resulting from ^-air interactions in the 
Earth's atmosphere is proportional to a v n- Fortuitously, it was recently showr^l 
that the flux of up-going charged leptons (UCL) per unit surface area produced by 
neutrino interactions below the Earth's surface is inversely proportional to o v n , as 
long as the neutrino absorption mean free path in Earth is small in comparison with 
the Earth's radius (i?©); i.e. for o v n ~ 2 x 10 _33 cm 2 (see Fig. EJ). This contrasts 
with the HAS rate proportional to a v ^ . Even when showering of the upgoing tau 
lcpton is included, the upgoing air-shower (UAS) rate still varies inversely with the 
cross-section in the range of interest, and can even exceed the HAS rate by several 
orders of magnitude, as displayed in Fig. |21 

Taken together, up-going and horizontal rates ensure a healthy total event rate, 
regardless of the value of o v n- Moreover, by comparing the HAS and UAS rates, 
the neutrino-nucleon cross section can be inferred at energies as high as 10 11 GeV or 
higher. This enables QCD studies at a minimum, and possibly discovery of a strong 
neutrino cross-section. o~ v n may also be determinable from a measurement of the 
angular distribution of UCL/UAS events, in addition to the approach comparing 
UAS and HAS rates. One expects the angular distribution of UCL to peak near 
cos6>p 0ak ~ (2 (p) R @ a vN )~ . 

Let us now examine the physics of upward showers in some detail. UHE neu- 
trinos are expected to arise from pion and subsequent muon decay. These flavors 
oscillate and eventually decohere during their Hubble-time journey. If |J7 T 3| — \U^\ 
are large, as inferred from the Super-Kamiokande data, then F Ut ss is expected. 
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Figure 2. Neutrino MFP displayed as a chord length in the earth. Cross-section values label the 
various MFP's. 



Because the energy-loss MFP for a t produced in rock or water is much longer than 
that of a muon or electron, the produced taus have a much higher probability to 
emerge from the Earth and to produce an atmospheric shower. Thus, the dominant 
primary for initiation of UAS "earth-skimming" events is the tau neutrino. 021 
Consider an incident tau neutrino whose trajectory cuts a chord of length I in the 
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Figure 3. The air shower probability per incident tau neutrino iJuAs/ F Ut tzA as a function of the 
neutrino cross section (eq.(7)). The incident neutrino energy is 10 20 eV and the assumed energy 
threshold for detection of UAS is E t ^ = 10 18 eV for curve 1 and 10 19 eV for curve 2. 



Earth. The probability for this neutrino to reach a distance x is P v {x) — e~ x / x ", 
where A^ 1 = <j v n P (the conversion from matter density to number density via 
iV^/gm is implicit). The probability to produce a tau lepton in the interval dx 
is j^-. The r produced at point x emerges from the surface with energy E T ~ 
E v e~( l ~ x )/ x T _ ^po p roc l UC e an observable shower, one requires P t ^uas = 9(A r + 
x — I), with A r = ((3 T Psr)^ 1 kL(E v /Eth)', Pt ~ 0.8 x 10~ 6 cm 2 /g is the exponential 
energy-attenuation coefficient and E t h is the minimum detectable energy. Taking 
the product of these conditional probabilities and integrating over the interaction 
site x we get the probability for a tau neutrino incident along a chord of length I 
to produce an UCL: 

P„ T ^(0 = f ^ = (e A ^ - 1) e~ 1 ^ . (5) 

The emerging tau decays in the atmosphere with probability 

Pd = 1 — cxp(— 2R^H/ct t I), where H 10 km parametrizes the height of the 
atmosphere. Thus, the probability for a tau-neutrino to produce an UAS is 

iV-UAs(0 = (1 - e' 2R ^ H /^ 1 ) P V ^ T {1) . (6) 

The fraction of neutrinos with chord lengths in the interval {1,1 + dl} is Trkrdl. 

Finally, including two further geometric factors, the solid angle 7r for a planar 
detector with hemispherical sky-coverage, and the tangential surface area A of the 
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detector, we arrive at the rate of UCL and UAS events: 

f 2R e> I dl 

^(uas, = F ^ A J q J^T Pu t ^t(VAS)(1) ■ (7) 

In Fig. 13 we show the number of expected UAS events per incoming neutrino as a 
function of the neutrino cross section. For comparison, we also show the number 
of expected HAS events per neutrino that crosses a 250 km held of view, up to an 
altitude of 15 km. It is clear that for the smaller values of the cross section, UAS 
events will outnumber HAS events, and vice versa. 

We give some examples of the UAS event rates expected from a smaller neutrino 
cross section at 10 20 eV, choosing a V N — 10 _33 cm 2 as an example. EUSO and OWL 
have shower-energy thresholds E t h ~ 10 19 eV corresponding to curve 2 in Fig. |3 
and apertures ~6x 10 4 km 2 and 3 x 10 5 km 2 , respectively. These detectors should 
observe F20 & n d 7Fzq UAS events per year, respectively (not including duty cycle); 
here F20 is the incident neutrino flux at and above 10 20 eV in units of km _2 sr _1 yr _1 , 
one-third of which are i/ T 's. Including showers from taus originating outside the 
field of view, and direct tau events, increases these rates, as does tilting space-based 
detectors toward the horizon to maximize the acceptance for events with smaller 
chord lengths and to allow more atmospheric path length for tau decay. The rates 
will also increase if E t h can be reduced. 

Recently, it was pointed out that production of the lightest supersymmctric 
particle (LSP) from annihilating dark matter may present a detectable flux of LSP- 
CRs.®A hig gsino-dominated LSP may have a cross-section as small as 10 2 times 
the SM neutrino cross-section. Such a tiny cross-section is a perfect example of 
what can be discovered and measured via either the UAS/HAS or the angular- 
dependence methods discussed here. 

4 Dispersion Relations: the High-Energy/Low-Energy Connection 

A neutrino primary is not subject to the GZK losses that limit cosmic nucleon 
propagation. However, the extrapolated neutrino-nucleon cross-section is expected 
to be be 10 -31 cm 2 , about 10 -6 too small to provide the neutrino with interactions 
high in our atmosphere. To explain the production of the observed EE cosmic ray 
events with neutrino primaries, some have postulated a new strong-interaction for 
neutrinos above ~ 10 19 eV. The idea of a strongly-interacting neutrino is not new, 
but recent developments in field theory and in gravity have given new motivation 
to such a picture. To mimic hadronically-induced air showers, the new neutrino 
cross section must be of hadronic strength, ~ 100 mb, above Eqzk = 5 x 10 19 eV. 
Simple perturbative calculations of single scalar or vector exchange cannot provide 
an acceptably fast growth of the cross-section with energy. El However, the modern 
thoughts on large TeV-scale cross-sections are much more imaginative. A plethora 
of new states, possibly growing exponentially in s or y/s, is motivated by precocious 
unification, low-scale string theory, and modes from additional space-dimensi ons 
accessible at y/a ~ TeV. Electroweak instantons are the most recent possibility. 

Direct limits on the EE neutrino cross-section are quite weak. The vertical 
column density of our atmosphere is X v = 1033 g/cm 2 . In terms of neutrino MFP 
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\ v , this may be written X v /\ u = ovn/1.6 mb. The horizontal slant depth 
is 36 times larger, leading to X^/X^ — a^jAA fib. Since penetrating events are 
not observed above us or to our side, the neutrinos must be interacting high in 
the atmosphere (large cross-section) or not interacting at all (small cross-section). 
Thus the cross-section range from ~ 20 fib to ~ 1 mb is excluded. More quantitative 
analyses give similar results. | 16 | 17 | 

A very interesting indirect limit on the EE neutrino cross-section is provided 
by dispersion relations. IAfl ' Dispersion relations are rigorous, nonperturbative, and 
model-independent. They limit the growth of the elastic neutrino amplitude at low 
energy due to any rising cross-section at higher energies. If new physics dominates 
the neutrino total cross-section with a value a* above the lab energy E*, then the 
dispersion relation determines the real part of the new strong-interaction elastic 
amplitude at lower energy E to be 'w G * ■ Remarkably significantly enhanced 
rates may occur for elastic vN scattering at an en ergy seven orders of magnitude 
lower than the onset of a new total cross-section. LL2l Such anomalous "low" en- 
ergy scattering may be observable in the neutrino beams available at Fermilab and 
CERN, and maybe in quasi-elastic e~p — > v e n scattering at HERA. 

How does this magic come about? Assuming only that the scattering amplitude 
is analytic, there results the following dispersion relation: El 

Re A ± {E) - Re A±(0) = -V f°° dE> ( a ^ E '^\ + ^(E',±) \ 

±y ' ±y ' 4tt J \E'(E'-E) E'(E' + E)J W 

where A±(E) are invariant v-N amplitudes, labeled by the nucleon helicity, and 
V denotes the principle value of the integral. Suppose the new physics dom- 
inates the neutrino-nucleon dispersion integral JBJ for E 1 > E* as hypothe- 
sized to explain the air showers observed above the GZK limit. Assuming that 
a* is independent of helicity and energy, and obeys the Pomeranchuk theorem: 
<j\o t (E, ±) — Otot (E, ±) o, the real part of the amplitude at energy E emerges: 

Rei±(_E)~Rci ± (0) + f ^ . (9) 

Z7T hi 

This result cannot be obtained in perturbation theory! Re ^4±(0) is nothing but the 
low energy limit of the weak interaction, ~ ^775- From this, we may immediately 
write down the ratio of the new amplitude to the SM amplitude: 



RcA(ff) ncw _ ( Ej\m GeV \ 

RcA(E)sm ~ \E*/10 18 eVJ { 100 mb 



(10) 



It is clear from (I10|l . and striking, that order 100% effects in the real elastic ampli- 
tudes begin to appear already at energies seven orders of magnitude below the full 
realization of the strong cross section. 

A promising observable consequence is available from the elastic cross section, 
obtained from the square of the elastic amplitude. The result ijHJ says that if the 
neutrino is strongly interacting at E* ~ 10 17 5 eV, then a factor of ten anomalous 
rise in the elastic cross-section is occurring at 100 GeV, a neutrino energy already 
available at Fermilab and CERN. Since the anomalous elastic cross-section grows 
quadratically with E, the anomalous event rate develops rapidly for E > 100 GeV. 
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Thus, the event sample of a future underground/ water/ice neutrino telescope opti- 
mized for TeV neutrinos could conceivably contain 1 000 times more elastic neutrino 
events than predicted by the SM; and a telescope optimized for PeV neutrinos may 
contain 10 9 more elastic events. Some of the wilder brane-world cross-sections 
proposed for the EE neutrino are ruled out by this dispersion result. 

There may be further tests of the strong-interaction hypothesis. If the neu- 
trino develops a strong-interaction at high energy, do not the electron and the 
other charged-lepton S[/(2)-doublet partners of the neutrinos also develop a similar 
strong-interaction? Is there new physics in the quasi-elastic e~p — > v e n scattering 
channel at HERA energies? A possible enhancement in the quasi-elastic channel 
cannot be deduced from dispersion relations. A separate calculation can be made, 
however, if certain aspects of the ne w hig h- energy strong-interaction are assumed. 
This is presently under investigation. 

5 Puzzles in the Extreme-Energy Cosmic Rays (EECRs) 

The discoveries by the AGASA, Fly's Eye, Haverah Park, and Yakutsk collabora- 
tions of air shower events with energies above the GZK cutoff challenge the SM of 
particle physics and the hot big-bang model of cosmology. Not only is the mecha- 
nism of particle acceleration for such EECRs controversial, but also the propagation 
of EECRs over cosmic distances is problematic. As has been mentioned, however, 
the HiRes experiment does not confirm the prior AGASA rate for events ~ 10 20 eV. 
The situation will remain murky until the Auger hybrid detector provides guidance. 

The famous Fly's Eye event occurred high in the atmosphere, whereas the event 
rate expected in the SM for early development of a neutrino-induced air shower 
is down from that of an electromagnetic or hadronic interaction by six orders of 

magnitude. On the other hand, Fig.0]presents the AGASA evidence that the arrival 

on 

directions of some of the highest-energy primaries are paired. Furthermore, a 
recent analysis of the arrival directions of the super-GZK events offers a tentative 
claim of a correlation with the directions of BL-Lac quasars. This correlation, 
if validated with future data, and the pairing data, argue for propagating cosmic 
particles which are charge neutral, stable, and have a negligible magnetic moment. 
The neutrino emerges as the only candidate among the known particles. 

Several non-neutrino solutions for the origin of the exceptional EECRs have been 
proposed. Magnetic caustics have been proposed as a focusing mini-lens bringing 
charged particles into pairs, and it is conceivable that some new physics at high 
energy could stabilize the neutron. With magnetic caustics one gets pairing but 
not pointing; with stabilized neutrons, one gets pairing and pointing. Generally, 
other proposed models are distinguishable from these and the neutrino scenario by 
the lack of pairing and pointing. 

The unexpected degree of small-scale clustering observed in the highest-energy 
AGASA cosmic ray day has motivated a derivation of analytic formulaJ22 which 
estimate the probability of random cluster configurations. The derived formulas 
offer a quick study of the strong potential of HiRes, Auger, and EUSO for deciding 
whether any observed clustering is meaningful or random. For detailed comparison 
to data, this analytical approach cannot compete with Monte Carlo simulations in- 
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Figure 4. AGASA sky-map of the arrival directions of the highest-energy cosmic rays (from the 
AGASA homepage). The five light-shaded disks contain doublets, and the one darker disk contains 
a triplet. Disk solid angles are n(9/2) 2 , with cone opening-angle 8 = 2.5°. 



eluding experimental systematics. Nevertheless, a recent Monte Carlo calculation^! 
of joint multiplet-probabilities shows reasonably good agreement with the analytic 
approach presented here. (The same MC analysis also cross-correlates different 
energy bins to suggest that pairing is emerging only above 10 19 eV.) The derived 
formulas here do offer two advantages over Monte Carlo techniques: (i) easy assess- 
ment of the significance of any observed clustering, and most importantly, (ii) an 
explicit dependence of cluster probabilities on the chosen angular bin-size. 

To derive the combinatoric formula for the probability of various event distri- 
butions in angle, imagine that the sky coverage consists of a solid angle £1 divided 
into N equal angular bins, each with solid angle ui ~ wO 2 steradian. Define each 
event distribution by specifying the partition of the n total events into a number 
mo of empty bins, a number mi of single hits, a number mi double hits, etc. The 
probability to obtain a given event topology is: 



P({ mj },n,N) 



N\ 



N'- 



(0!) m ° (l!) m i (2!) m 2 (3!) m 3 



(11) 



The variables in the probability are not all independent. The partitioning of 
events is related to the total number of events by X)j=i 3 x m i = n i an( l 1° the 
total number of bins by Ylj=o m i = ^ ■ Because of these constraints, one infers 
that the process is not described by a simple multinomial or Poisson probability 
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distribution. It is useful to use these constraints to rewrite the probability (|1 If) as 

p(w,n,iv)= (i2) 



where 



W-iv(-) -. (13) 

In the n <C iV limit, is expected to approximate the mean number of j-plets, 
and eq. I|12|l becomes roughly Poissonian. As an approximate mean, rnj defined 
in eq. I)13|l provides a simple estimate of cluster probabilities due to chance for the 
7i<C N case. 

Two large- number limits of interest are N ^S> n ^> 1, and n > N 3> 1. With bin 
numbers typically ~ 10 3 , the first limit applies to the AGASA, HiRes, and Auger 
experiments; the second limit becomes relevant for the EUSO experiment. Valid 
when N 3> n 3> 1, one has: 



P({mJ,n, N) 



n 



(14) 



where r = (AT — mo)/n « 1, and the prefactor (nearly unity) is 

P = e^ (l ^" ll) . (15) 



The non-Poisson nature of Eq. (|14fl is reflected in the factorials and powers of r in 
the exponents, and the deviation of the prefactor from unity. 

In the case where n > N 3> 1, higher j-plets are common and the distribution 
of clusters can be rather broad in j. We may write rnj in the approximate form: 

_ / ATS / en Y + 5 , . 

Extremizing this expression with respect to j, one learns that the most populated 
j-plct occurs near j ~ n/N. Combining this result with the broad distribution 
expected for large n/N, one expects clusters with j up to several x to be common 
in the EUSO experiment. 

Shown in Fig. is an assessmenJ23l of the AGASA-plets, five doublets plus a 
triplet, obtained using formula l|14|l . Two features of the figure are noteworthy. The 
first is the rather extreme sensitivity of the statistical significance to the angular 
binning size. AGASA claims a 2.5° resolution, which puts the significance of their 
clusters at 10 -3 . If the resolution were 3° (2.0°), the significance would be a factor 
of six weaker (fifteen stronger) . The second feature is the factor of a few error made 
when Poisson statistics are blindly applied. 

The explicit dependence of random clustering probabilities on angular bin-size 
presented here may prove quite useful in the future. If clustered events originated 
from a common source and traveled without bending, then the experimental angular 
resolution is the optimal bin-size. On the other hand, if clustering results from 
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Figure 5. Exact (solid) and Poissonian (dashed) inclusive probabilities (from ref. 24) for five 
doublets and one triplet in the 58-event AGASA sample. 



magnetic focusing, then the angular size of magnetic caustics may be the relevant 
bin-size. If clustering results from density fluctuations in the Galactic halo, then 
the angular size of these fluctuations may be the optimal bin-size. Since photons 
are not bent by magnetic fields whereas protons are bent, the optimal bin-size for 
photon-initiated events is likely smaller than that for proton-initiated events. The 
analytic formulas are easily applied to any chosen angular bin-size. 

6 Z-bursts 

A rather conservative and economical solution proposed to solve the super-GZK 
mysteries is the Z-burst mechanism. Here, EECR neutrinos scattering resonantly 
on the cosmic neutrino background (CNB) predicted by Standard Cosmology, to 
produce Z-bosons. ESI These Z-bosons in turn decay to produce a highly boosted 
"Z-burst" , containing on average twenty photons and two nucleons above Eqzk (see 
Fig. EJ. The photons and nucleons from Z-bursts produced within 50 to 100 Mpc 
of earth can reach earth with enough energy to initiate the air-showers observed at 
~ 10 20 eV. The energy of the neutrino annihilating at the peak of the Z-pole is 
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V 

COSMIC RAY 




Figure 6. Schematic diagram showing the production of a Z-burst resulting from the resonant 
annihilation of a cosmic ray neutrino on a relic (anti)ncutrino. If the Z-burst occurs within the 
GZK zone (~ 50 to 100 Mpc) and is directed toward the earth, then photons and nucleons with 
energy above the GZK cutoff may arrive at earth and initiate super-GZK air-showers 



E % = 2^ = HeV/mi)ZeY - 

The resonant-energy width is narrow, reflecting the narrow width of the Z-boson: 
at FWHM AE R /E R ~ T z /M z = 3%. The mean energies of the - 2 baryons 
and <~ 20 photons produced in the Z decay are easily estimated, when the Z-burst 
energy is averaged over the mean multiplicity of 30 secondaries in Z-decay. The 
photon energy is reduced by an additional factor of 2 to account for their origin in 
two-body 7T° decay. 

In the simplest approximation, the spectrum of arriving nucleons is 

dN 1 dN dD ^ , 



weiler'YITP: submitted to World Scientific on February 2, 2008 16 



from sources uniformly distributed out to 

In ( NE GZK \ 

D GZK ~X^=±, (19) 

with a pileup at Eqzk resulting from all primaries originating beyond this dis- 
tance. The E~ x spectrum extends from Eqzk out to the maximum nucleon energy 
~ E R /30 ~ 10 21 ( °^ cV ) eV. More realistic simulations including energy- loss pro- 
cesses, cosmic expansion, and boosted Z-boson fragmentation functions give a softer 
spectrum, but a characteristic feature of the Z-burst mechanism remains that the 
super-GZK spectrum is considerably harder than the sub-GZK spectrum having 
power law index -2.7. 

The necessary conditions for the viability of this model are then, a neutrino 
mass scale of the order 0.1 to 1 eV. and a sufficient flux of neutrinos at > 10 21 
e V. The second condition seems challenging, while the first is quite natural in 
view of the recent neutrino oscillation data. F its t o atmospheric neutrino data yield 
m v > V^™atm ~ 0.05 eV. A recent estimate^! of the total neutrino mass in the 
Universe, based on the distribution of large-scale structures, is Y^j m j ~ 2 eV. Most 
recently, the WMAP collaboration ambitiously combined their new CMB data with 
matter-distribution spectra to deduce'23 J^j m j SS> 0.71 eV. However, omission of 
the suspect Lyman-alpha data from the analysis returns one to Ylj m j ^ 1 e V, 
with dependence on priori^. And so it appears that the neutrino mass is squeezed 
to lie within just the 0.1 to 1.0 eV range most beneficial to the Z-burst model! 

What is not known is whether Nature has provided the large neutrino flux at 
energy E R to allow an appreciable event rate in future EECR detectors. It is 
conceivable, although unlikely, that the flux is so large that present EECR events 
are initiated by Z-bursts. A recent analysiJ^of this possibility gave a best fit with 
m, y = 0.26^0 ^4 eV, nicely consistent with the WMAP bound. Another analysiJ^l 
fits the EECR spectrum down to the ankle with Z-burst generated events and a 
neutrino mass of y/ <5m 2 tm ~ 0.07 eV, again in accord with the WMAP bound. 
The flux requirements for the Z-burst mechanism can be ameliorated if there is an 
overdensity of relic neutrinos, as would happen if (i) there was a significant chemical 
potential, or (ii) neutrinos were massive enough to cluster in "local" structures 
such as the Galactic SuperCluster. Large chemical potentials have been ruled out 
recentljl^, and this exclusion is confirmed by the WMAP data. Local clustering 
has been studiedp^ with the conclusion being that a significant overdensity on the 
SuperCluster scale requires a neutrino mass in excess of 0.3 eV. Such a mass is 
marginally allowed by the new WMAP limit. 

The large neutrino flux required at E R to explain the AGASA super-GZK event 

rate with the Z-burst mechanism has engendered much debate. While the direct 

JT7I 

neutrino flux limits do not preclude such a large flux, limits on the associated 
diffuse gamma-ray flux disfavor most source mechanisms. A wide range of possi- 
ble sources have been critically reviewed recently. I^4 | 35 | g ven t ne so-called hidden 
sources (which emit only neutrinos) seem problematical. EH On the other hand, 
new po ssib ilities have emerged for production of neutrino fluxes. Wake-fields in 
plasmaJ23 may accelerate particles to energies of 10 23 eV, and a special class of 
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blazars may beam neutrinos at us. ^Hl "Mirror matter" topological d efects may 
decay to mirror neutrinos which then may oscillate into active neutrino J^^ I Re- 
combination of the strong magnetic fields surrounding black holes offers another 
acceleration mechanism that is still in the infancy of its exploration. EDI 

7 Summary and Acknowledgments 

The next ten years will exploit neutrino eyes to see the extreme Universe. The- 
ory and experiment are rapidly converging to usher in the neutrino window to 
astronomy, and the astronomy window into neutrino physics. Fruitful discoveries, 
probably beyond anything anticipated today, await us. 
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